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. Abstract. A simple model of nucleons coupled to angular momentum zero (s-pairs) occupying the 

valance shell of a semi-magic nuclei is considered. The model has a separable, orbit dependent 
pairing interaction which dominates over the kinetic term. It is shown that such an interaction leads 
to an exactly solvable model whose (0 + ) eigenstates and energies can be computed very easily with 
the help of the algebraic Bethe ansatz method. It is also shown that the model has a supersymmetry 
which connects the spectra of some semimagic nuclei. The results obtained from this model for 
the semimagic Ni isotopes from 58 Ni to 68 Ni are given. In addition, a new and easier technique for 
calculating the energy eigenvalues from the Bethe ansatz equations is also presented. 
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We consider the s-pairs of the same kind of nucleons (either protons or neutrons) oc- 
cupying the valance shell of a semi-magic nuclei and interacting with an orbit dependent 
separable pairing force as described by the Hamiltonian 
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The orbit dependence of the pairing force is described by the dimensionless constants cj 
while the overall strength of the pairing term against the kinetic term is measured by the 
constant \G\ which has the dimension of energy. In writing the Hamiltonian (Q~|), we use 
the quasispin operators which are given by |Q]] 

Sj = (Sj)'= Z (- 1 ) y "" fl J» fl }-« = * E (^Uj m + a]_ m aj- m -l). (2) 

m>0 Z m>0 

Quasispin operators obey the angular momentum commutation relations 

[Sf,Sj f ] = 28 jf $ [tf,Sp] = ±8 jf Sf, (3) 

i.e., one has an angular momentum algebra for each orbit j such that those angular 
momenta corresponding to different orbits commute with one another. Note that Eq. © 
implies —Q.j/2 < S® < Q.j/2 where Clj = j + 1/2 is the maximum number of s-pairs 
which can occupy the orbit j. This tells us that the quasi-spin algebra corresponding to 
the level j is realized in the Q.j/2 representation. 

Here, we are interested in the limit of the pairing Hamiltonian in which all the 
single particle energy levels become degenerate. In this degenerate limit, the pairing 



Hamiltonian has two interesting properties which are discussed in this talk: 1) It is 
exactly solvable and 2) it is invariant under a supersymmetry transformation which 
connects the spectra of different isotopes or isotones. Note that the pairing Hamiltonian 
is also exactly solvable away from the degenerate limit (i.e., with different single particle 
energy levels) but only in the presence of two orbits [2y . From this solution, we know 
that the supersymmetry mentioned above is broken away from the degenerate limit. 

We begin by noting that in the degenerate limit, the one body term of the pairing 
Hamiltonian is reduced a constant proportional to the total number of pairs which can 
be discarded. Therefore, in the degenerate limit, the pairing Hamiltonian becomes 

H D = -\G\C^cp+)(^c f S f ). (4) 

J f 

This tells us that the degenerate limit of the pairing Hamiltonian can also be viewed as 
the strong pairing limit. 

Exact eigenvalues and eigenstates of the Hamiltonian given in Eq. © were obtained 
in Refs. O, 0, (see Ref. [@] for the eigenvalues of the quantum invariants of this 
Hamiltonian). Following these references, we introduce the operators 

5+ W = I 1 ? i2 s t and s~(*) = E i r 12 S J- (5) 

J J J J 

These operators respectively create and annihilate a single pair of nucleons in the valance 
shell in a certain distribution over the single particle energy levels described by the 
complex parameter x. The eigenstates of the pairing Hamiltonian in the degenerate limit 
can be written in terms of these operators. For example, in Ref. [7] Talmi has showed 
that, in the presence of one pair of nucleons in the valance shell, 

5+(0)|0) E = -\G\ 1 £a j \c j \ 2 (6) 

j 

is an eigenstate-eigenvalue pair of the Hamiltonian ©. Additional eigenstates with one 
pair of nucleons can be found by considering a more general Bethe ansatz state in the 
form of 5 + (x)|0). In fact, using the quasi-spin angular momentum algebra given in Eq. 
©, it is easy to show that 

S + (x)\0) E = (7) 

is an eigenstate-eigenvalue pair of the Hamiltonian © if the parameter x is a solution of 
the Bethe ansatz equation 
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Note that Eq. ([8]) may have several distinct solutions in which case we have a zero energy 
eigenstate in the form of Eq. © for each solution. It is worth to emphasize at this point 
that all the energies are calculated with respect to the core and that we have ignored the 
one body term in the Hamiltonian. 



These considerations generalize to the case of several pairs occupying the valance 
shell. Let us denote the maximum number of pairs which can occupy the valance shell 
by Nmax- If the shell is occupied by N pairs such that N < N max /2 then the pairing 
Hamiltonian given in Eq. © has some eigenstates with nonzero energy and some 
eigenstates with zero energy, similar to the one pair case discussed above 0,0,0]. The 
nonzero energy eigenstates together with their energies are given by 

^(0)^( Zl )...^( ?Y _ 1 )|0) E = -\G\ fanjlcf- N f-) (9) 

V j fe=i Zk J 

where the parameters zi,Z2, ■ . . , zn-1 are solutions of the set of Bethe ansatz equations 
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for m= l,2,...,N— 1. Similarly, the zero energy eigenstates are given by 

S + (x l )S + (x 2 )...S + (x N )\0) E = (11) 
where the parameters x\,X2, ■ ■ ■ ,xjv are solutions of the Bethe ansatz equations 
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for m = 1, 2, . . . , N. Note that the Bethe ansatz equations given above for both the non- 
zero energy eigenstates and the zero energy eigenstates may have several distinct sets 
of solutions. In such a case, each solution set gives us an eigenstate-eigenvalue pair 
when substituted in the relevant formulas. In particular, if the zero energy Bethe ansatz 
equations have more than one solutions then the zero energy level is degenerate. 

Let us now consider a valance shell which is more than half full, i.e., N > N max /2. 
One can, in principle, extend the formulas (l9~lfT2~l) for N > N max /2. But in this case one 
frequently encounters one of the following two problems: 1) The Bethe ansatz equations 
given above yield no solutions or 2) the Bethe ansatz equations do yield solutions but 
the corresponding states vanish identically. For this reason, it is more reasonable to use 
hole pairs instead of nucleon pairs to describe a valance shell which is more than half 
full. 

The pairing Hamiltonian given in Eq. © is invariant under a supersymmetry trans- 
formation which transforms the spectrum of a nuclei with Af pairs of valance nucleons 
(N < Nmax/2) to the spectrum of a nuclei with N — 1 pairs of valance holes. In order to 
describe this symmetry, let us first introduce the particle-hole transformation operator 



f = exp l- ix y^—J-\ (13) 




transforms the empty valance shell |0) into the fully occupied valance shell |0) and 
the pair creation- annihilation operators into each other, i.e., 

f t|0) = |0) f f Sff = Sj f f S°jf = -S°j (14) 



Therefore, it transforms a state containing N particle-pairs into states containing N hole- 
pairs. However, there is no particle-hole symmetry in the problem described by the 
pairing Hamiltonian as can be easily verified by showing that [Hp, 7 1 ] ^ 0. On the other 
hand, if we define the operators 



S + = (£c*S+)f B- = f\^cjSJ 



(15) 



then it is easy to show that B + and B commute with one another 

[£+,£]=() 

and that the pairing Hamiltonian given in Eq. (HJ) can be written as 

H D = - GB B + = -\G\B + B~. 



(16) 
(17) 
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This tells us that the problem described 

by the pairing Hamiltonian in the de- i c g- 
generate limit possesses a quantum me- 
chanical supersymmetry. Since the op- 
erators B + and B commute with each 
other, the Hamiltonian Ho is equal to 
its own supersymmetric partner. As a 
result, if the state is an eigenstate 
of the Hamiltonian Hd with a nonzero 
energy, then the state B | l//) is also an 

eigenstate of H D with the same energy. N par tj C | e . pa j rs 

Zero energy states, however, are annihi- 
lated by B and do not transform under 
the supersymmetry. It is easy to see from 

the definition of B given in Eq. £[5]) FIGURE 1. Supersymmetry transformation, 
that if the state 1 l//) has /V particle-pairs, 

then the state B \ ty/) has N —I hole-pairs. Because the operator B first annihilates one 
particle-pair and then turns the remaining N —I particle-pairs into hole-pairs. Consider, 
for example, the nonzero energy eigenstate with one pair of particles given in Eq. ©. 
Using Eqs. ([141) and (fTTT) . one can easily show that 



N-1 hole-pairs 



zr£ + (o)|o)<x |0) 



(18) 



Consequently, the states S + (0) |0) and |0) (i.e., the fully occupied valance shell) have the 
same energy eigenvalue given by Eq. ©. On the other hand, the zero energy eigenstates 
described by Eqs. © and © are annihilated by B . Similarly, the nonzero energy 
eigenstates with /V pairs of nucleons described in Eqs. © and (flOl) transform as 



B-S + (0)S + (zi)...S + (zn-i)\0)ocS-(zi)...S-(zn-i)\0) 



(19) 



As a result, the states S+(0)S+(zi) . . .S+(z N -\)\0) and S~(zi) ...§'(z N -i)\0) have 
the same energy given in Eq. ©. Since /V < N max /2, the former state corre- 
sponds to a nucleus which has less than half full or half full valance shell 



whereas the later state corresponds to a nucleus which has more than half 
full valance shell (See Figure [[])■ Note that the zero energy states with N 
pairs of nucleons described by Eqs. ([111) and (PT21) are annihilated by B . 

As an application of this formalism, 

c 1/2 = 0.31 5 2p 1/2 i e t us consider the even-even Ni iso- 

■^i/2 = 1 topes from 58 Ni to 68 Ni for which the 
neutron and proton shells are closed at 
N = Z = 2S (i.e., the core is the dou- 

° 5 ' 2 - 5 ", bly magic 56 Ni nucleus). The extra neu- 

5/2 * trons in 58 Ni- 68 Ni occupy the single par- 
0777 9 n ticle energy levels shown in Figure El 
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q _2 The occupation probability amplitudes 
shown in this figure are taken from Ref. 
Jil] and the value of the pairing strength 
\G\ is chosen so as to obtain the ground 
state energy of 58 Ni (with respect to 56 Ni 
FIGURE 2. The valance shell for 58 Ni- 68 Ni. core) in agreement with the experimen- 
tal data. Figure[3] shows the exact energy 
eigenvalues of the pairing Hamiltonian given in Eq. © calculated with the formalism 
outlined above. The zero energy states of 58 Ni and 60 Ni isotopes are doubly degenerate 
as indicated by double lines in the figure. All the energies are given in MeV and are 
relative to the Ni core. We use the supersymmetry to obtain the spectra of the nuclei 
with more than half full shell (i.e., 64 Ni- 68 Ni). For example, the spectrum of Ni is the 
same as the spectrum of 60 Ni except that the zero energy states are absent for 64 Ni. See 
Ref. [4] for a comparison of these results with the experimentally available values. 

With the Bethe ansatz formalism outlined above, one can easily calculate the exact 
numerical values of the pairing energies. However, to obtain closed analytical expres- 
sions for the pairing energies is relatively difficult in this formalism. Because this would 
require the triple effort of solving the Bethe ansatz equations analytically, substituting 
the solutions into the energy expression given in Eq. © and then making the necessary 
simplifications. In the last part, we present an easier method of calculating the energy 
eigenvalues. Although this method is currently valid for only two orbits, it can possibly 
be extended to include more than two orbits. 

In the presence of only two orbits, the problem of calculating the nonzero energy 
eigenvalues can be greatly simplified by using a method proposed by Stieltjes [9]. Using 
this method, one can transform the problem of solving the equations © and finding the 
energy eigenvalues for /V nucleon-pairs into the problem of finding N th order polynomial 
solutions of the following differential equation fl 1 QTI : 

p » (z) -(% + %lV (z)+ (ggjZgO PNiz) = (20) 
V z z—i ) z{z- 1) {z — c) 



where a = -N, /3 = N — Qj ] — Clj 2 - 1 and 
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FIGURE 3. The exact energy eigenvalues of the pairing Hamiltonian for Ni- Ni in MeV relative to 
the 56 Ni core. The degeneracies of the zero energy states of 58 Ni- 60 Ni are indicated by the double lines. 



Note that a similar but slightly different method was used in Ref. ^\ in order to obtain 
solutions of the zero energy Bethe ansatz equations. Eq. (|20l) is known as Heun's 
differential equation and it admits polynomial solutions only for certain values of the 
parameter q. One can use the Frobenius method (i.e., try a series solution) and then find 
the conditions on q such that the series will terminate at the N th order. This gives us the 
allowed values of the energy eigenvalues for N pairs of nucleons because q and E are 
related by Eq. (f2TT) . This method significantly reduces the amount of calculation needed 
to compute the energy eigenvalues of the pairing Hamiltonian given in Eq. © II 1011 . 
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